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24 ABSTRACT 

25 Positive-sense (+) RNA viruses encode RNA-dependent RNA polymerases (RdRps) 

26 essential for genomic replication. With the exception of the large nidoviruses, such as 

27 coronaviruses (CoVs), RNA viruses lack proofreading and thus are dependent on RdRps to 

28 control nucleotide selectivity and fidelity. CoVs encode a proofreading exonuclease in 

29 nonstructural protein 14 (nspl4-ExoN), which confers a greater-than 10-fold increase in fidelity 

30 compared to other RNA viruses. It is unknown to what extent the CoV polymerase (nspl2- 

31 RdRp) participates in replication fidelity. We sought to determine whether homology modeling 

32 could identify putative determinants of nucleotide selectivity and fidelity in CoV RdRps. We 

33 modeled the CoV murine hepatitis virus (MHV) nspl2-RdRp structure and superimposed it with 

34 solved picornaviral RdRp structures. Fidelity-altering mutations previously identified in 

35 coxsackie virus B3 (CVB3) were mapped onto the nspl2-RdRp model structure and then 

36 engineered into the MHV genome with [nspl4-ExoN(+)] or without [nspl4-ExoN(-)] ExoN 

37 activity. Using this method we identified two mutations conferring resistance to the mutagen 5- 

38 tluorouracil (5-FU): nspl2-M61 IF and nspl2-V553I. For nspl2-V553I we also demonstrate 

39 resistance to the mutagen 5-azacytidine (5-AZC) and decreased accumulation of mutations. 

40 Resistance to 5-FU, and decreased number of genomic mutations, was effectively masked by 

41 nspl4-ExoN proofreading activity. These results indicate that nspl2-RdRp likely functions in 

42 fidelity regulation and that, despite low sequence conservation, some determinants of RdRp 

43 nucleotide selectivity are conserved across RNA viruses. These results also indicate that, with 

44 regards to nucleotide selectivity, nspl4-ExoN is epistatic to nspl2-RdRp, consistent with its 

45 proposed role in a multi-protein replicase/proofreading complex. 

46 
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47 IMPORTANCE 

48 RNA viruses have evolutionarily fine-tuned replication fidelity to balance requirements 

49 for genetic stability and diversity. Responsibility for replication fidelity in RNA viruses has been 

50 attributed to the RNA-dependent RNA polymerases, with mutations in RdRps for multiple RNA 

51 viruses shown to alter fidelity and attenuate virus replication and virulence. Coronavimses 

52 (CoVs) are the only known RNA viruses to encode a proofreading exonuclease (nsp 14-ExoN), as 

53 well as other replicase proteins involved in regulation of fidelity. This report shows that the CoV 

54 RdRp (nsp 12) likely functions in replication fidelity, that residue determinants of CoV RdRp 

55 nucleotide selectivity map to similar structural regions of other unrelated RNA viral 

56 polymerases, and that for CoVs, the proofreading activity of the nsp 14-ExoN is epistatic to the 

57 function of the RdRp in fidelity. 
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58 INTRODUCTION 

59 RNA virus replication results in the incorporation of a relatively high number of 

60 mutations, ranging from 10" 4 to 10" 6 mutations per site per round of replication (1-5). It is thought 

61 that low-fidelity replication is largely responsible for the capacity of RNA viruses to evolve 

62 rapidly and adapt to new host species and ever changing environmental pressures (6-8). The 

63 RNA-dependent RNA polymerase (RdRp) is central to the replication of RNA viruses and is a 

64 key regulator of nucleotide selectivity and fidelity (9, 10). Recent studies of CVB3, polio, HIV-1 

65 and other viruses demonstrated that viable viruses are only recoverable within a four-fold range 

66 of RdRp fidelity (11-14). In most cases altered RdRp fidelity decreases fitness relative to WT 

67 viruses; this has been demonstrated for changes as small as a 1.2 fold difference in the 

68 accumulation of mutations (12, 14-16). Despite having as low as no amino acid identity outside 

69 of conserved motifs (11-14, 17-19), all described polymerase structures (including RdRps) 

70 resemble a “cupped right hand” with fingers, palm and thumb domains (20). The fingers form a 

71 channel that allows entry of the template RNA and ribonucleotide triphosphates (rNTPs) and 

72 assist in proper positioning of incoming nucleosites in the active site (21). The palm contains the 

73 active site, and the thumb functions in contacting exiting nascent RNA (21-23). However, there 

74 is diversity in the viral proteins that encode for RdRps; additional domains are often present that 

75 encode a variety of functions such as methyltransferase, endonuclease, polyribonucleotidyl- 

76 transferase, guanylyltransferase, membrane targeting, protein-protein binding or protein-RNA 

77 binding activities (24-26). 

78 Coronaviruses (CoVs) infect a wide array of species and have emerged as highly 

79 pathogenic human pathogens twice in this century, first with Severe Acute Respiratory 

80 Syndrome coronavirus (SARS-CoV) in 2003 (27) and second with Middle East Respiratory 
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81 Syndrome coronavirus (MERS-CoV) in 2012 (28). CoVs, and other large nidoviruses, replicate 

82 with higher fidelity than all other known (+)-RNA viruses (29, 30). CoVs also encode for the 

83 largest known RNA virus genomes ranging from 27-34 kb (31, 32), and increased fidelity in 

84 CoVs is likely required for maintenance of these large genomes (14). CoV genomes encode 16 

85 non-structural proteins (nspl-16), several of which are known or predicted to function in fidelity 

86 regulation, including nspl4-ExoN, a 3’-5’ exoribonuclease, and nsplO, a modulator of nspl4- 

87 ExoN activity (33, 34). Mutating the DE of the nspl4-ExoN active site to AA inactivates the 

88 exoribonuclease, yielding nspl4-ExoN(-) viruses, and nspl4-ExoN(-) viruses exhibit a greater- 

89 than 10-fold increase in mutation frequency (29, 35-37). Recent evidence demonstrates that 

90 nspl4 directly interacts with the CoV RdRp encoded in nspl2 (nspl2-RdRp) (38), but the effect 

91 of this interaction on nucleotide selectivity and overall fidelity regulation is not known. There are 

92 no solved structures for any CoV nspl2-RdRp, but the presence of conserved RdRp motifs and 

93 modeling of the C-terminal half of nspl2 predicts an RdRp domain that is structurally similar to 

94 other RNA viruses (39, 40). 

95 The demonstrated function of nsp 14-ExoN in high fidelity CoV replication raises the 

96 question of whether and how nspl2-RdRp participates in fidelity regulation. We sought to 

97 determine whether nspl2-RdRp can modulate nucleotide selectivity independently or in 

98 association with the proofreading nsp 14-exonuclease. We modeled the RdRp domain of CoV 

99 nsp 12 on coxsackievirus vims B3 (CVB3) and poliovirus polymerase structures and predicted 

100 residues important for fidelity based on prior result from those virus systems. Substitution 

101 mutations at these residues were introduced in the isogenic recombinant genome of the |3-CoV, 

102 murine hepatitis vims (MHV-A59). We demonstrate that two of these mutations, nspl2-V553I 

103 and nspl2-M61 IF, confer resistance to the mutagen 5-fluorouracil (5-FU) and one, nspl2- 
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104 V553I, also results in resistance to the mutagen 5-AZC and demonstrates a decreased 

105 accumulation of mutations. Increased mutagen resistance and decreased accumulation of 

106 mutations was only observed in viruses with an inactivated ExoN, demonstrating that nsp 14- 

107 ExoN proofreading activity is epistatic to the nucleotide selectivity of nspl2-RdRp. In this paper 

108 we define epistasis as a situation where the phenotype of one gene or viral protein masks the 

109 phenotype of genetic variants of another viral protein. This result is consistent with a primary 

110 role for nsp 14-ExoN in error recognition and removal. However, introduction of RdRp mutations 

111 within the WT-MHV background decreased fitness relative to WT. Together the results suggest 

112 that nspl2-RdRp shares common determinants of nucleotide selectivity with RdRps from other 

113 RNA virus families. Further, the CoV RdRp has likely evolved to function in cooperation with 

114 nsp 14-ExoN rather than independently. 
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115 MATERIALS AND METHODS 

116 Virus and cell culture. Murine delayed brain tumor (DBT) cells (41) and baby hamster kidney 

117 21 cells expressing the murine hepatitis virus (MHV) receptor (BHK-R) (42) were maintained at 

118 37°C in Dulbecco’s modified Eagle medium (DMEM; Gibco) supplemented with 10% FBS 

119 (Invitrogen), penicillin, streptomycin (Gibco) and amphotericin B (Corning). BHK-R cells were 

120 further supplemented with 0.8mg/mL of G418 (Mediatech). All virus work was performed using 

121 recombinant WT-MHV strain MHV-A59 (GenBank accession number AY910861 (42)). 

122 

123 Sequence analysis and homology modeling of CoV MHV nspl2-RdRp. The MHV RdRp 

124 domain structure was generated using Phyre2 online program (43) using nspl2 residues 385-887 

125 that correspond to the reported SARS-CoV nspl2-RdRp model (40). The structural model was 

126 compared to the X-ray crystal structures of coxsackivirus B3 (CVB3) and poliovirus (PDB 

127 accession numbers 3DDK and 1RA7 respectively) using the Pymol Molecular Graphics System 

128 (Schrodinger, LLC). ClustalX multiple-sequence alignments were generated using the program 

129 Mac Vector. 

130 

131 Cloning, recovery and verification of mutant viruses. Quick-change mutagenesis was used to 

132 generate point mutations in individual MHV genome cDNA fragment plasmids using the 

133 previously described MHV infectious clone reverse genetics system (42). Mutant viruses were 

134 recovered in co-cultured BHK-R and DBT cells following electroporation of in vitro transcribed 

135 genome RNA in BHK-R cells. All viruses that included nspl4-ExoN(-) mutations were 

136 generated using the F fragment previously described (35). Before use in viral recovery all 

137 mutagenized plasmids were fully sequenced (GenHunter Corporation, Nashville, TN) to assure 
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138 no additional mutations were introduced. We also sequence verified engineered mutations in 

139 recovered viruses. Viruses in the nspl4-ExoN(-) background took between 84 and 96 hours to 

140 reach around 80 percent involvement in syncytia for a P0 stock in contrast to viruses in the WT 

141 background which were frozen at 24-48 h p.i. PI working stocks were made by infecting D9s at 

142 an MOI 0.01 and freezing when 80 percent involved in syncytia, approximately 24 h p.i for WT 

143 viruses and 36 h p.i for nspl4-ExoN(-) viruses (2-3 rounds of replication). 

144 

145 Compounds and drug sensitivity studies. 5-fluorouracil (5-FU) was obtained from Sigma, and 

146 prepared as 200mM stock solutions in DMSO. 5-azacytidine (5-AZC) was also obtained from 

147 Sigma and prepared as 50mM stocks solutions in water. Sub-confluent DBT cells were pre- 

148 treated for 30min with DMEM with the indicated concentrations of 5-FU or DMSO, 5-AZC or 

149 media alone. Treatment was removed, and inocula added and allowed to adsorb for lhr at 37°C. 

150 Inocula was then removed and media with drug or DMSO was returned. Infection proceeded for 

151 24hrs for WT or 32hrs for nspl4-ExoN(-) viruses, when supernatants were acquired, frozen and 

152 titered by plaque assay as previously described (35). 

153 

154 Virus replication and RNA synthesis assays. Sub-confluent DBT cell monolayers in triplicate 

155 were infected at an MOI of 0.01 PFU/cell. Virus was allowed to adsorb for 30min when inocula 

156 were removed and the cells were washed 2X with PBS, followed by addition of pre-warmed 

157 media. For replication kinetics assays samples were taken at various time points post infection. 

158 Titering was performed by plaque assay as previously described (35). For analysis of RNA 

159 synthesis, total infected-cell RNA was obtained by TRIzol reagent (Invitrogen) at various times 

160 post infection and two-step RT-qPCR was performed as previously described (44). 
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161 

162 Determination of Specific Infectivity. Sub-confluent DBT cells were pre-treated for 30min 

163 with DMEM with indicated concentrations of 5-FU or DMSO alone. Treatment was removed, 

164 and inocula added and allowed to adsorb for lhr at 37°C. Inocula was then removed and media 

165 with drug or DMSO was returned. Infection proceeded for 20hrs for WT or 24hrs for nspl4- 

166 ExoN(-) viruses then supernatants were acquired, frozen and titered by plaque assay as 

167 previously described (35). Supernatants were also used for RNA genome isolation by adding 

168 lOOuL supernatant to 900uL TRIzol reagent (Invitrogen), chloroform extracting by phase 

169 separation and using the aqueous layer in the PureLink Mini RNA kit (Ambion) as per the 

170 manufacturers protocol. One-step RT-qPCR was performed as below and the ratio of PFU to 

171 genomes of the supernatant was determined. 

172 

173 One-step RT-qPCR for determining supernatant genome copies for specific infectivity 

174 assay. An RNA standard was generated using MHV A fragment (42) to generate a 931 

175 nucleotide RNA. First cDNA was generated by PCR using the primers: forward 5’- 

176 T A AT ACGACT C ACT AT AGGGGGCT AT GT GG ATT GTTGT GG-3 ’ which begins with a T7 

177 promoter and reverse 5’-AATTCTTGACAAGCTCAGGC-3’. RNA for the standard curve was 

178 then generated using mMessage mMachine T7 kit (Ambion). An agarose gel with 1% bleach was 

179 run and a ~900nt band was observed. RNA was purified using RNeasy Mini kit (Qiagen). 

180 Dilutions of the standard curve were made from 10 3 to 10 s genome equivalents for use in assay 

181 as needed. Primers and probes for one-step RT-qPCR were purchased from BioSearch Tech. 

182 Probe is 5’ FAM labeled and 3’ BHQ-1 labeled with the sequence 5’- 

183 TTCTGACAACGGCTACACCCAACG-3 ’ and made up to 5uM in nuclease free water. The 
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184 primers used were forward 5’-AGAAGGTTACTGGCAACTG-3’ and reverse 5’- 

185 TGTCCACGGCTAAATCAAAC-3’. Reactions were set-up on ice with enzyme added last. 

186 Final volume for reactions was 20uL with 150nM probe, 900nM each primer, 2uL sample RNA 

187 and lOuL 2X ToughMix, one-step, low ROX enzyme mix (Quantas) used per reaction. Samples 

188 were plated in duplicate and run on the Applied Biosciences 7500 Real-Time PCR System with 

189 the conditions 55C for lOmin, 95C for 5min, 95C for 30sec and 60C for lmin with the last two 

190 steps repeated 40X. The standard curve was graphed and genomes per mL determined. 

191 

192 Competitive fitness of mutant viruses. Sub-confluent DBT monolayers were co-infected at a 

193 total MOI of 0.01 PFU/mL with RdRp mutant viruses in the nspl4-ExoN(-) background and 

194 nspl4-ExoN(-) at either a 1:1, 1:9 or 9:1 ratio. When 50-70% of the monolayer was involved in 

195 syncytia, total RNA was harvested. RNA was then reverse transcribed using SuperScriptHI 

196 (Invitrogen) per the manufactures protocol and amplicons were generated using primers designed 

197 to cover the region including the codons for both the V553 and M611 residues. Amplicons were 

198 sent for sequencing and electopherograms were analyzed using MacVector. 

199 

200 Passage reversion analysis. Triplicate monolayers of sub-confluent DBT cells were infected 

201 with an initial MOI of 0.01 PFU/mL of nspl2-V553I and nspl2-M61 IF viruses in both the WT 

202 and nspl4-ExoN(-) background. Viruses were then blind passaged in triplicate for 5 passages. 

203 Total RNA was sequenced across a 1.7kb region of nspl2-RdRp that included both nspl2-RdRp 

204 mutations. Electropherograms were analyzed using MacVector. 

205 

206 Preparation of amplicons for deep sequencing of full viral genomes. Sub-confluent DBT 
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207 cells were infected at an MOI of 0.01 PFU/mL with nspl2-V553I or nspl2-M61 IF in either the 

208 WT or nspl4-ExoN(-) background, nspl4-ExoN(-) alone or WT alone. Infections were allowed 

209 to progress for 20 hours then RNA was isolated. RNA was reverse transcribed using 

210 SuperScriptlll (Invitrogen) per the manufactures protocol, and 12 amplicons were generated to 

211 cover the whole genome and processed as described previously (44). 

212 

213 Deep sequencing sample preparation and analysis. Amplicons were subsequently purified via 

214 a nucleospin PCR purification kit (Macherey-Nagel), quantified by picogreen, fragmented 

215 (Fragmentase) and prepared using the lllumina NextSeq500 Mid Output 150 cycle kits following 

216 the standard protocols. Sequences were obtained with an lllumina NextSeq500 machine. 

217 Sequencing runs were analyzed using the previously published ViVan bioinformatics pipeline 

218 (45). Briefly, the pipeline performs quality filtering and adaptor cleaning was done using fastq- 

219 clipper (http://hannonlab.cshl.edu/fastx_toolkit/index.html). The 150-nt reads were aligned to the 

220 reference sequence with a maximum 2 mismatches per read, using BWA (46) and processed 

221 using SAMTOOLS (47) to obtain the nucleotide/base calling at each position. The ViVAn 

222 pipeline then identifies statistically significant variants above the background noise due to 

223 sequencing error, calculated for each nucleotide site: for each position throughout the viral 

224 genome, base identity and their quality scores are gathered. Each variant allele's rate is initially 

225 modified according to its covering read qualities based on a maximum likelihood estimation, and 

226 tested for significance using a generalized likelihood-ratio test. Additionally, an allele confidence 

227 interval is calculated for each allele. In order to correct for multiple testing, Benjamini-Hochberg 

228 false-discovery rate of 5% was set. In all experiments, a minimum coverage of 3000X reads was 

229 obtained and the background error frequency at every nucleotide site was always below 0.0001. 
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230 For analysis we use a conservative frequency cut off of 0.01 consistent with previous studies (48- 

231 50). 

232 

233 Statistical analysis. Statistics were applied as described in figure legends using GraphPad Prism 

234 6 software (La Jolla, CA). The number of replicates performed for each experiment is similarly 

235 listed in each respective figure legend. Finally, some of the data was normalized to controls; the 

236 GraphPad Prism 6 software also performed this analysis. 
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237 RESULTS 

238 Homology modeling of MHV nspl2-RdRp polymerase core domain predicts putative 

239 fidelity determinants. Mutations that alter nucleotide selectivity have been identified across 

240 multiple RNA virus RdRps (3, 11, 51-53); however, whether these residues are conserved across 

241 virus families is unknown. We sought to determine whether residues within nspl2-RdRp that are 

242 structurally homologous to known RNA vims fidelity determinants would have similar effects on 

243 nucleotide selectivity when introduced into the MHV background. To do this, we modeled the 

244 structure of MHV nspl2-RdRp using Phyre2 software (43). A series of nspl2-RdRp truncations 

245 was assessed, and the highest confidence model was used for further study. This region 

246 corresponded with a published model for the SARS nspl2-RdRp (40) and included residues 385- 

247 887 of the MHV nspl2 protein, referred to here as the RdRp core domain (Fig. la,b). Deletion of 

248 the CoV specific domain (residues 1-384), and a small C-terminal portion of the thumb domain 

249 (888-928), was required to establish this high-confidence model (Fig. la). The model was 

250 resolved by highest probability similarity to human rhinovirus serotype lb (PDB ID: lXR5a), 

251 rabbit hemorrhagic disease vims (PDB ID: 1KHV) and EV71 (PDB ID: 3N6M). The Phyre2 

252 confidence, i.e. the probability of true homology, for the RdRp core domain to these structures 

253 was >99% while the percent identity was only 14-20%. Having generated a structural model for 

254 the MHV nspl2-RdRp core domain we next sought to predict residues involved in nucleotide 

255 selectivity. The nspl2-RdRp core domain model was aligned with the solved structure of 

256 coxsackie virus B3 (CVB3, PDB ID: 3DDK) using PyMol (Fig. lc). A series of CVB3 RdRp 

257 mutations have been shown to result in decreased fidelity (10, 11, 54). The CVB3 fidelity 

258 determinants were compared with the MHV nspl2-RdRp core domain model. Those that aligned 

259 well structurally and by amino acid similarity were further investigated: MHV nspl2-V553, 
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260 M611, W613, A621, Y649 and K794 (Fig. 1). Finally, the nspl2 aa sequences of 27 different 

261 a-, (3- and y-CoVs were aligned, including SARS-CoV and MERS-CoV. All six identified 

262 residues were conserved across these CoVs (Fig. Id). Analysis of similarity and types of residues 

263 in the picornaviruses were then used to determine the specific amino acid changes that would be 

264 introduced at the identified MF1V residues. The resulting substitution mutations were engineered 

265 in the isogenic cloned MHV genome: nspl2-V553A/I, M611F, W613Y, A621G, Y694F1/W and 

266 K794R (Table 1). 

267 

268 Recovery of mutant viruses in the MHV nspl4-ExoN(+) and nspl4-ExoN(-) isogenic 

269 backgrounds. We next tested whether viable viruses could be recovered with substitutions at 

270 the identified residues. Virus recovery was attempted a maximum of 3 times, resulting in 

271 recovery of 6 of the 8 mutant viruses in the WT background: nspl2-V553I, M61 IF, W613Y, 

272 A621G, Y649H, and K794R. The time required for recovery of mutant viruses in the WT 

273 background ranged from 24 to 48 hours. No other mutations were identified across nspl2 

274 sequences in recovered viruses. The A621G mutant was not further studied as it demonstrated 

275 rapid primary reversion even in the recovery (P0) supernatant. Since our goal was to understand 

276 the relationship of nspl2-RdRp and nsp 14-ExoN in fidelity regulation we additionally attempted 

277 recovery of the WT background viable mutants in the setting of inactivated nsp 14-ExoN (nsp14- 

278 ExoN(-)). In contrast to mutant viruses recovered in the WT background, we only recovered 2 of 

279 the 5 mutants in the nspl4-ExoN(-) background: nspl2-V553I/nspl4-ExoN(-) and nspl2- 

280 M61 lF/nspl4-ExoN(-) (Table 1). The time to recovery for nspl2-V553I/nspl4-ExoN(-) was 84 

281 hours and for nspl2-M61 lF/nspl4-ExoN(-) was 96 hours. Working stocks of all viruses were 

282 made by infecting D9s with an MOI of 0.01 and recovering stocks at around 24 h p.i. for WT 
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283 viruses, or between 32 and 48 h p.i. for mutant viruses in the nspl4-ExoN(-) background. 

284 Therefore 4-10 replication cycles were required to generate stocks. Working stocks were 

285 sequenced to verify that the introduced mutations were still present. 

286 

287 Resistance of recovered mutant viruses to the base analog 5-fluorouracil. We next tested our 

288 panel of recovered mutant viruses for resistance to the RNA mutagen 5-fluorouracil (5-FU). 5- 

289 FU has been used with picornaviruses, influenza viruses, vesicular stomatitis viruses and others 

290 to reflect changes in fidelity based on increased or decreased sensitivity to incorporation and 

291 virus inhibition (11, 30, 49, 55). WT CoVs (MF1V, SARS-CoV) are resistant to 5-FU, while 

292 nspl4-ExoN(-) mutants are profoundly sensitive to 5-FU inhibition, consistent with nspl4- 

293 ExoN-mediated removal of misincorporated 5-FU. The effect of 5-FU on DBT cell viability was 

294 previously tested with no effect observed up to 400uM (30). We compared WT-MF1V with the 

295 nspl2 mutants in both the WT and nspl4-ExoN(-) backgrounds. There was no significant change 

296 in sensitivity to 5-FU compared to WT for any of the nspl2-RdRp mutant viruses in the WT 

297 background at up to 120uM 5-FU although the mutation Y649FI does appear to decrease 

298 resistance slightly (Fig. 2a). In contrast, the nspl2-M61 lF/nspl4-ExoN(-) and nspl2- 

299 V553I/nspl4-ExoN(-) mutant viruses both were significantly less sensitive to 5-FU than nspl4- 

300 ExoN(-) alone, with both populations persisting when treated with 120uM 5-FU where nspl4- 

301 ExoN(-) was not detectable beyond 80uM 5-FU (Fig. 2b). These data demonstrate that both 

302 nspl2-RdRp mutations V553I and M61 IF confer resistance to 5-FU. This suggests a couple 

303 possibilities: that it is not possible to increase the exclusion of 5-FU beyond the high level 

304 dictated by nspl4-ExoN or that selectivity for native nucleotides over 5-FU is in fact increased 

305 by nspl2 mutations, but at a low level that is not detectable by changes in virus titer. 
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306 Three of the mutations (nspl2-Y649H, W613Y, and K794R) were viable in the nspl4-ExoN(+) 

307 background, but failed to grow in the absence of proofreading. The observation that they only 

308 grew in an nspl4-ExoN(+) background indicated that the mutations retained sufficient 

309 polymerase function to support virus replication, but it was critically dependent on having 

310 proofreading functionality, suggesting these mutations may have given rise to low fidelity 

311 variants. Unfortunately, the low titer from the nspl4-ExoN(-) background precluded direct 

312 sequencing analysis and we therefore cannot definitively show this is the case. 

313 

314 Replication kinetics of nspl2-V553I and nspl2-M611F mutant viruses in the WT and 

315 nspl4-ExoN(-) backgrounds. Since we were interested in mutations that potentially confer 

316 altered fidelity, we prioritized the nspl2-V553I and nspl2-M61 IF mutant viruses for further 

317 analysis. We next sought to determine how the nspl2-V553I and nspl2-M61 IF viruses 

318 replicated in comparison to their isogenic background (Fig. 3). In the wild-type background both 

319 mutant viruses had slightly delayed exponential replication, but eventually reached similar peak 

320 titers comparable to WT. In contrast, in the nspl4-ExoN(-) background both nspl2-V553I and 

321 nspl2-M61 IF mutant viruses displayed similar replication kinetics to the isogenic nspl4-ExoN(- 

322 ) background. We also assessed RNA synthesis for nspl2-V553I and nspl2-M61 IF by RT- 

323 qPCR. Measured genomic RNA levels were consistent with the virus replication kinetics data, 

324 and we observed delayed and decreased genome RNA synthesis in the WT background after 

325 multiple rounds of replication (Fig. 3E). However, at early time points there was no difference in 

326 RNA accumulation suggesting that decreased RNA is a result of steps post RNA synthesis (Fig. 

327 3C). In the nspl4-ExoN(-) background, RNA synthesis levels were indistinguishable from the 

328 nspl4-ExoN(-) background for both nspl2-V553I and nspl2-M611F, with no additional RNA 
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329 synthesis defects detectable (Fig. 3D and F). The results, along with the ability to recover several 

330 of the nspl2-RdRp mutants only in the nspl4-ExoN(+) background, support the hypothesis that 

331 nspl4-ExoN and nspl2-RdRp may have an epistatic relationship. However, for replication 

332 kinetics the effects of nspl2-V553I and nspl2-M61 IF are not observable in the presence of an 

333 inactive nspl4-ExoN, demonstrating that the replication phenotype of an inactive npsl4-ExoN is 

334 epistatic to replication variants encoded in nspl2-RdRp. 

335 

336 Specific infectivity of nspl2-V553I and nspl2-M611F. Having identified two mutant viruses 

337 with resistance to 5-FU we wanted to further test whether resistance was due to decreased 

338 incorporation of the mutagen. Measurement of specific infectivity has been useful for 

339 determining lethal mutagenesis for MHV and other RNA viruses (30, 56). We tested the nspl2- 

340 M61 IF and nspl2-V553I mutants, in both the WT and nspl4-ExoN(-) backgrounds, for changes 

341 in specific infectivity when infected at an MOI of 0.01 and treated with 5-FU (Fig. 4). Both 

342 nspl2-M61 IF and nspl2-V553I resulted in an increased ratio of infectious particles (PFU/ml) to 

343 total particles (RNA genomes) in the nspl4-ExoN(-) background (Fig. 4b). Similarly, in the WT 

344 background nspl2-M611F demonstrated an increase in the ratio of infectious particles to RNA 

345 genomes (Fig. 4a). Thus specific infectivity may be a more sensitive measure of lower level 

346 changes in nucleotide selectivity in the setting of nspl4-ExoN(-). 

347 

348 Fitness cost of nspl2-V553I and nspl2-M611F in WT and nspl4-ExoN(-) backgrounds. In 

349 multiple RNA viruses, including CoVs, both increased and decreased fidelity have been reported 

350 to have a fitness cost (14). We therefore sought to determine whether nspl2-V553I or nspl2- 

351 M61 IF conferred any cost in fitness, defined as the ability to directly compete during co- 
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352 infection. In the WT background, both nspl2-V553I and nspl2-M61 IF demonstrated delays in 

353 replication and impaired RNA accumulation. However, since there were not any observed 

354 additional defects in replication or RNA synthesis for nspl2-V553I and nspl2-M61 IF when 

355 introduced in the nspl4-ExoN(-) background (Fig. 3b, d), we tested for any additional fitness 

356 cost of nspl2-V553I/nspl4-ExoN(-) or nspl2-M61 lF/nspl4-ExoN(-) viruses compared with 

357 nspl4-ExoN(-) alone. When co-infected with nspl4-ExoN(-) at ratios from 1:9 to 9:1, nspl2- 

358 V553I/nspl4-ExoN(-) maintained the input ratio compared with nspl4-ExoN(-) (Fig. 5a). A 

359 small advantage to nspl2-V553I/nspl4-ExoN(-) was observed when the co-infected cultures 

360 were treated with 60uM 5-FU, consistent with a conferred advantage for 5-FU resistance (Fig. 

361 5a). Thus there appeared to be no additional fitness cost of nspl2-V553I/nspl4-ExoN(-) 

362 compared to nspl4-ExoN(-) alone. In contrast, nspl2-M61 lF/nspl4-ExoN(-) was not able to 

363 compete with nspl4-ExoN(-) at any ratio (Fig. 5b). Treatment with 60uM 5-FU again favored 

364 nspl2-M61 lF/nspl4-ExoN(-); however, even then the percent of the population made up of 

365 nspl2-M61 lF/nspl4-ExoN(-) only remained at around 30% of the population when initially 

366 given a 9-fold advantage (Fig. 5b). 

367 We next tested whether the relative differences in fitness cost resulted in selective pressure for 

368 reversion of nspl2-V553I and nspl2-M611F (Fig. 6). DBT cells were infected with nspl2- 

369 V553I, nspl2-V553I/nspl4-ExoN(-), nspl2-M611F, or nspl2-M61 lF/nspl4-ExoN(-) at an 

370 initial MOI of 0.01. After 5 passages the viruses were analyzed for retention of original 

371 mutations using di-deoxy (Sanger) sequencing. The nspl2-V553I mutation was stable after 

372 passage in both WT and nspl4-ExoN(-) backgrounds, maintaining the mutated AUU codon (Fig. 

373 6a), suggesting minimal selective pressure on that nucleotide, codon, or amino acid. In contrast, 

374 the nspl2-M61 IF mutation demonstrated significant change over passage in both the WT or the 
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375 nspl4-ExoN(-) backgrounds. The original mutation, UUC, was no longer the majority codon in 

376 the nspl4-ExoN(-) background and was less than 52% of the population in all WT background 

377 lineages (Fig. 6b). The nspl2-M61 lF/nspl4-ExoN(-) population resulted in a mix of two 

378 nucleotide changes resulting in reversion to AUG (methionine) (<68% of the population), single 

379 nt changes that result in mutation it to Leu (<55% of the population), a somewhat smaller but 

380 still hydrophobic residue, or retaining a Phe substitution (<27% of the population). Thus the 

381 fitness cost of nspl2-M61 IF results in significant selective pressure for changes at that residue 

382 during passage in absence of 5-FU. 

383 

384 Resistance of recovered mutant viruses to the base analog 5-azacytidine. Having shown that 

385 both nspl2-M61 IF and nspl2-V553I mutations conferred resistance to 5-FU in the nspl4- 

386 ExoN(-) background but that nspl2-M611F reverted quickly and was at a fitness disadvantage in 

387 both the WT and nspl4-ExoN(-) background, we next wanted to test whether the mutations 

388 conferred resistance specifically to 5-FU or broadly to base analogs and therefore were likely 

389 determinants of fidelity. A study by Arias et al demonstrated that for foot-and-mouth disease 

390 virus (FMDV) resistance to a specific mutagen can result from a point mutation while conferring 

391 the opposite overall fidelity (57). We therefore tested for resistance to the additional mutagen 5- 

392 azacytidine (5-AZC). Similar to the 5-FU results, the nsp 12-V553I/nsp 14-FxoN(-) mutant virus 

393 was more resistant to 5-AZC than nspl4-ExoN(-) alone, maintaining approximately one log 

394 higher titers from 20-50uM 5-AZC (Fig. 7). However, the nspl2-M61 lF/nspl4-ExoN(-) mutant 

395 virus showed no difference in resistance to 5-AZC compared with nspl4-ExoN(-) alone (Fig. 7). 

396 These data suggest that nspl2-V553I is likely a fidelity determinant where as nspl2-M61 IF 
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397 confers specific resistance to 5-FU and has an unknown overall fidelity. This result is also 

398 consistent with the rapid reversion observed for the nspl2-M611F mutation. 

399 

400 The nspl2-V533I mutation results in a decrease in the accumulation of mutations. Flaving 

401 shown that the nspl2-V553I and nspl2-M61 IF mutations resulted in resistance to 5-FU, and that 

402 nspl2-V553I additionally conferred resistance to 5-AZC, we sought to directly determine 

403 whether either of these mutations resulted in a change in the number of mutation accumulated in 

404 viral RNA. DBT cells were infected with wild type, nspl2-V553I, nspl2-M611F, nspl4-ExoN(- 

405 ), nspl2-V553I/nspl4-ExoN(-) or nspl2-M61 lF/nspl4-ExoN(-) at an MOI of 0.01 and RNA 

406 was collected at 20 h. pi. These samples were then prepared for Illumina next generation 

407 sequencing (NGS) across the full genome and analyzed using the ViVan analysis pipeline (45). 

408 Mutations present as 1% or more of the population were graphed by frequency and position in 

409 the genome. Engineered mutations are depicted with colored dots. Most non-engineered new 

410 mutations were present at 10% or less of the population and were distributed across the genome 

411 with no detectable “hot-spots” (Fig. 8a-f). For nspl2-V553I in the WT background, no difference 

412 was observed in the number of mutations accumulated to 1% or greater of the population 

413 compared to WT alone. In contrast, in the nspl4-ExoN(-) background, nspl2-V553I was 

414 associated with a 1.7-fold decrease in the frequency of mutations compared to nspl4-ExoN(-) 

415 alone (Fig. 8g), again with no change in the distribution of mutations across the genome (Fig. 8c, 

416 d). These results are consistent with both the 5-FU and 5-AZC data in suggesting increased 

417 fidelity. Results from the nspl2-M61 IF mutant viruses were more complicated. To our surprise, 

418 the nspl2-M61 lF/nspl4-ExoN(-) virus fully reverted at both engineered nspl2-M61 IF 

419 nucleotides during the low MOI infection, resulting in a viral population that was nspl4-ExoN(-) 
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420 alone (Fig. 8f). This reversion made the results for the nspl2-M61 lF/nspl4-ExoN(-) virus 

421 interpretable. Flowever, since accumulation of mutations to over 1% of the population is a 

422 combination of all replication cycles from initial recovery to the final sample (roughly 11 total 

423 for nspl2-M61 lF/nspl4-ExoN(-)) we included this data in Figure 8g and 8h. In the WT 

424 background the nspl2-M61 IF mutations were still present as 100% of the population and 

425 resulted in a 1.93 fold increase in the total number of accumulated mutations. We observed only 

426 a slight increase in the number of mutations accumulated in the nspl2-M61 lF/nspl4-ExoN(-) 

427 sample over those of nspl4-ExoN(-), which was not surprising seeing as the M61 IF mutation 

428 was no longer present. Neither sample appeared to have mutations concentrated in specific 

429 locations across the genome suggesting the accumulation of mutations was due to random 

430 generation of mutations rather than strong selection in particular locations or proteins (Fig. 8e 

431 and f). Of note, one mutation in nsp3 of the nspl4-ExoN(-) sample and a mutation in nsp3, nspl3 

432 and the E proteins of nspl2-V553I/nspl4-ExoN(-) reached nearly 100% of the sample 

433 population. None of these mutation were present in the fragments used for recovery and the nsp3 

434 mutation which arose in the nspl4-ExoN(-) population was not present prior to the final low 

435 MOI infection. It is possible that these mutations provide some benefit to these viruses since they 

436 were fixed so rapidly in the population. We next determined whether either of the mutations 

437 resulted in a change in the types of mutations occurring during replication (Fig. 8h). Consistent 

438 with our previous studies, there were differences in the types of mutations incorporated when 

439 comparing WT and nspl4-ExoN(-) backgrounds (30). However, the addition of nspl2-V553I or 

440 nspl2-M61 IF did not alter these patterns in either WT or nspl4-ExoN(-) backgrounds. Thus, 

441 nspl2-V553I results in an overall decrease in the accumulation of mutations over passage while 

442 nspl2-M61 IF seems to increase the number of mutations accumulated. These results confirm by 
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443 sequence analysis the results from the 5-FU and 5-AZC resistance experiments for nspl2-V553I, 

444 specifically that the effects of nspl2-V553I are dependent on inactivation of nspl4-ExoN for 

445 their detection and that nspl2-V553I likely confers broad resistance to the incorporation of 

446 incorrect nucleotides. In contrast, these results further increase the complexity of the nspl2- 

447 M61 IF mutation in relation to incorporation of nucleotides and their analogs. We conclude that 

448 nspl2-M61 IF confers resistance to 5-FU but that this resistance is not likely to be due to a broad 

449 resistance to the incorporation of alternate nucleotides. 
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450 DISCUSSION 

451 RdRp structures of divergent RNA vimses are structurally conserved and likely have 

452 common determinants of activity in the finger, palm and thumb domains. Positive strand RNA 

453 virus polymerases appear to utilize a common palm-domain based mechanism for active site 

454 closure (23) and associated molecular determinants of fidelity in different RdRp domains have 

455 been proposed based on biochemical and mutagenesis studies (10). To date there are no solved 

456 crystal structures of any CoV RdRp and thus direct comparison with other virus RdRp structures 

457 has not been possible. Further, regulation of CoV fidelity is likely dependent on multiple 

458 proteins, including the RdRp and proofreading ExoN. Thus it was not clear that a CoV would 

459 phenotypically exhibit effects from mutating fidelity-determining residues located in the RdRp 

460 itself. In this study we sought to determine whether we could use structure and mutagenesis data 

461 from distantly related RNA viruses to identify determinants of CoV nspl2-RdRp fidelity. Our 

462 results suggest that CoV RdRps do in fact participate in fidelity regulation, at residues 

463 orthologous to those in the picornaviruses, and likely in other RNA virus RdRps. The results also 

464 define, for the first time, a CoV RdRp determinant that increases resistance to multiple 

465 mutagens, decreases the accumulation of mutations over time and so, likely increases overall 

466 fidelity. Additionally, the data in this paper suggest that CoV RdRp-mediated increased fidelity 

467 is only detectable when nsp 14-ExoN is inactive, and is only partially compensating for the loss 

468 of nsp 14-ExoN high fidelity. Finally, both the nspl2-V553I and M61 IF mutations confer a 

469 replication cost in the WT background but only nspl2-M611F confers a fitness disadvantage in 

470 the nspl4-ExoN(-) background. Together the results suggest that nspl4-ExoN proofreading 

471 activity is epistatic to nspl2-RdRp fidelity but that in contrast the replication defects of an 

472 inactive nspl4-ExoN are epistatic to replication defects in nspl2-RdRp. 
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473 

474 Determinants of nucleotide selectivity and fidelity in Co Vs may be conserved with other 

475 RNA viruses. Picornavinis functional RdRps contain only RdRp domains (58, 59). This is not 

476 the case for many viral RdRps, CoVs included (26, 39). In addition to the predicted RdRp core 

477 domain, all CoV nspl2 proteins contain a “CoV-specific” domain of over 350 amino acids at the 

478 N terminus of the protein. A nucleotidyltransferase activity was identified in this CoV specific 

479 domain and has been shown to be important in SARS-CoV replication, but the specific function 

480 in replication remains to be determined (39). The N-terminal CoV specific domain could not be 

481 modeled due to lack of evolutionary homologues of known structure but, the predicted RdRp 

482 core domain could be modeled with high confidence using bioinformatic approaches. The 

483 structures which provided the best models for the CoV RdRp were from picornaviruses, 

484 including enterovirus 71, foot and mouth disease virus and coxsackie virus B3. This allowed for 

485 direct alignment and comparison of the known fidelity determinants in CVB3 with the MHV 

486 nspl2 core domain, across both fingers and palm domains. The recovery of mutations in these 

487 structurally conserved residues and their participation in CoV nucleotide selectivity supports the 

488 hypothesis that there are determinants of base specificity conserved between Co Vs and distantly 

489 related RNA viruses, specifically the picornaviruses. This supports the idea that all RdRps 

490 function similarly due to their structural conservation and despite the low level of sequence 

491 similarity and attached domains (54). It also suggests that the CoV RdRp domain folds in a 

492 manner similar to other RdRps, likely separate from the CoV-specific domain. 

493 RdRp fidelity and nucleotide selectivity has been investigated extensively in 

494 picornaviruses, especially poliovirus (12, 15, 51, 53, 60, 61). However, even between 

495 picornaviruses, the impact on fidelity of changes at identical or similar residues can vary 
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496 dramatically (10). Mutations at the same residues in the RdRp ofpoliovims and CVB3 affect 

497 fidelity differently; poliovirus mutations generally result in increased fidelity and CVB3 

498 mutations in decreased fidelity (11). Specifically, the structurally orthologous residue to nsp 12- 

499 M611 in poliovirus is F230 and in CVB3 it is 1230; when F230 is mutated to an lie in poliovirus 

500 it results in an increase in fidelity, but when 1230 in CVB3 is mutated to a Phe, Trp or Val it 

501 results in a decrease in fidelity (10). Similarly, the orthologous residue to nspl2-V553 in CVB3 

502 is 1176; when 1176 is mutated to a Val in CVB3 it results in a decrease in fidelity (10, 11). 

503 Although the specific substitutions differ, when nspl2-V553I and nspl2-M61 IF are introduced 

504 into the MHV genome the resulting 5-FU resistance, and for nspl2-V553I a likely increase in 

505 fidelity, mimics results seen in poliovirus rather than CVB3. Flowever, nspl2-M61 IF seems to 

506 accumulate more mutations overtime than controls and many of the predicted mutations were 

507 only recoverable in the WT background, therefore changes in nucleotide selectivity and fidelity 

508 were not tested, so it remains possible that the mutations that were non-viable in the nspl4- 

509 ExoN(-) background had decreased fidelity and that nspl4-ExoN(-) treated with 5-FU defines 

510 the error threshold for CoVs. 

511 

512 Coronavirus nspl2-RdRp and nspl4-ExoN cooperate to optimize both fidelity and 

513 replication kinetics. In addition to nspl2-RdRp, CoVs encode nsp 14-ExoN, which functions as 

514 a proofreading enzyme (29, 30, 33, 35). Beyond nsp 14-ExoN, there are additional CoV-encoded 

515 nsps that potentially contribute to overall fidelity, such as: the small molecule modulator of 

516 nspl4-ExoN encoded in nsplO (33, 62), nsp7 and nsp8 that together function as an elongation 

517 factor (63) and a primase (64, 65), and nspl3 that functions as a helicase (66). These proteins 

518 interact with each other and likely function as a multi-protein replication/fidelity complex (34, 
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519 38). The results of this study show that while determinants in nspl2-RdRp are likely capable of 

520 increasing fidelity, those changes are only detectable in the setting of loss of nsp 14-ExoN 

521 proofreading, and do not completely compensate for the impaired fidelity associated with nsp 14- 

522 ExoN(-). If altering nucleotide selectivity determinants within nspl2-RdRp no longer 

523 significantly affect the overall nucleotide selectivity of the WT virus, then the evolutionary 

524 pressure on nspl2-RdRp may be more heavily weighted toward other aspects of replication, such 

525 as speed. Recent evidence suggests that RdRp fidelity and speed have an inverse relationship 

526 (10, 67). So, if fidelity regulation by nspl4-ExoN occurs more rapidly than correct nucleotide 

527 selectivity by the RdRp then CoV nspl2-RdRp may have been selected specifically for 

528 replication speed. Our data may support this hypothesis, as in the WT background both nspl2- 

529 V553I and nspl2-M61 IF resulted in increased replication lag phases (Fig. 3A and E), though 

530 early RNA synthesis was not observably different in our system (Fig. 3C). Flowever, this was not 

531 seen in the nspl4-ExoN(-) background which could be explained by our data indicating that 

532 replication is already slowed when nsp 14-ExoN is inactivated. One possible mechanism for this 

533 could be that the nspl4-ExoN(-) protein is trying to remove incorrect nucleotides and stalling 

534 replication due to its inability to do so. In this case nspl4-ExoN(-) would be epistatic to nspl2- 

535 RdRp in relation to speed, thereby obscuring decreases in replication speed caused by mutations 

536 in nspl2-RdRp itself. 

537 

538 Conclusion. Our results support the hypothesis that determinants of nucleotide selectivity are 

539 conserved across viral orders, identify the first likely increased fidelity determinant for CoV 

540 nspl2-RdRp, and demonstrate that nspl4-ExoN proofreading activity is epistatic to nspl2-RdRp 

541 nucleotide selectivity. Knowing that some fidelity determinants may be conserved across viral 
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542 orders is an exciting discovery, as many fidelity determinates identified so far have resulted in 

543 attenuation (11, 15, 54, 68-70). It would be interesting to determine whether nspl2-V553I is also 

544 attenuated. We would predict that it would be, based on the fitness cost in vitro. However, there 

545 is also no clear increase in nucleotide selectivity or fidelity in the WT background. This may be 

546 due to a very minimal change that is just not measureable even with the deep sequencing 

547 technology used in this study. Alternatively, it could be that any change is simply overwhelmed 

548 by the high fidelity of intact ExoN proofreading. In any case, these mutations or other changes at 

549 these residues may allow for selection of viruses that replicate with normal kinetics in vitro and 

550 in vivo, yet confer attenuation in an animal setting. We know that the nspl4-ExoN(-) mutations 

551 confer genotypically and phenotypically stable attenuation in vivo (36). However, the concept of 

552 a high level mutator as a mechanism for attenuation in live viruses may be problematic. The 

553 identification of increased fidelity mutations in the RdRp that can partially, or potentially 

554 completely, compensate for the fidelity impairment of nspl4-ExoN(-) viruses, may allow for 

555 development of approaches that can benefit from the stability of the nspl4-ExoN(-) mutator 

556 phenotype while allowing more stability to the input genomes. Finally, these results combined 

557 with those from previous work (33, 44) suggest that CoVs encode at least three proteins involved 

558 in fidelity (nspl2-RdRp, nspl4-ExoN and nsplO) supporting the assembly of a multi-protein 

559 replicase / fidelity complex as described previously (38). This increases the importance of 

560 establishing a biochemical model of the multi-protein complex to directly test the interactions of 

561 fidelity determinants as well as potential inhibitors of each or all of these functions. 
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760 FIGURE LEGENDS 

761 Figure 1. Homology modeling of CoV nspl2-RdRp and identification of residues that 

762 potentially regulate fidelity based on CVB3 structure. Phyre2 software was used to model A) 

763 a subsection of the MHV nspl2-RdRp core domain (expanded from nspl2-RdRp full schematic). 

764 This B) modeled MHV RdRp structure was aligned with the C) solved CVB3 RdRp structure. 

765 Residues chosen for site-directed mutagenesis (A and B) were selected by comparing previously 

766 determined fidelity altering mutations of picornavirus RdRps. Amino acid alignments D) across 

767 CoVs show all residues are completely conserved. 

768 

769 Figure 2. Resistance of MHV nspl2-RdRp mutant viruses to 5-fluorouracil in the WT and 

770 nspl4-ExoN(-) background. Domain location of mutations is indicated by: fingers (blues) and 

771 palm (reds). DBT cells were pre-treated with different concentrations of 5-FU for 30 min. 

772 Treatment was removed and cells were infected with indicated viruses in the A) WT background 

773 or B) nspl4-ExoN(-) background at an MOI of 0.01. Media containing 5-FU was replaced 30 

774 min p.i. Vims samples were taken at 24 (WT) or 32 (nspl4-ExoN(-)) h p.i. and titer was 

775 determined by plaque assay. Data represents 3 independent experiments, each with 2 replicates. 

776 Error bars represent SEM (*P < 0.05 by Wilcoxon test). 

777 

778 Figure 3. Replication kinetics of MHV nspl2-RdRp mutant viruses. Mutation location is 

779 indicated by: fingers (blues) and palm (red). DBT cells were infected with the viruses indicated 

780 in the A, C and E) WT background or B, D and F) nspl4-ExoN(-) background at an MOI of 0.01 

781 PFU/cell. Supernatant aliquots were taken at indicated times p.i. and titer determined by plaque 
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782 assay. Total RNA was taken at indicated times p.i. and RT-qPCR was performed. Data 

783 represents 3 independent experiments. Error bars represent SEM. 

784 

785 Figure 4. Specific infectivity is increased in both nspl2-V553I and nspl2-M611F mutants. 

786 DBT cells were pre-treated with increasing concentrations of 5-FU for 30 min. Treatment was 

787 removed and cells were infected with indicated viruses in the A) WT background or B) nspl4- 

788 ExoN(-) background at an MOI of 0.01. Media containing 5-FU was replaced 60 min p.i. Virus 

789 samples were taken at 20 and 24 h p.i. respectively. Titer was determined by plaque assay and 

790 number of supernatant genomes were determined using one-step RT-qPCR. Data represents 2 

791 independent experiments, each with 3 replicates. Error bars represent SEM (*P < 0.05 by 2way 

792 ANOVA using the Bonferroni correction for multiple comparisons). 

793 

794 Figure 5. Competitive fitness analysis in the WT and nspl4-ExoN(-) Background. DBT cells 

795 were pre-treated with media alone or containing 60uM or 300uM 5-FU for 30 min. Treatment 

796 was removed and cells were co-infected, at total MOI of 0.01, with WT and A) nspl2-V553I or 

797 B) nspl2-M611F at a ratio of 1:9 or nspl4-ExoN(-) and C) nspl2-V553I/nspl4-ExoN(-) or D) 

798 nspl2-M61 lF/nspl4-ExoN(-) at a ratio of 9:1, 1:1 or 1:9. Media alone or containing 60uM or 

799 300uM 5-FU was replaced 30 min p.i. Total RNA was taken at 20 h p.i. for WT competition (A 

800 and B) or 24 h p.i. for nspl4-ExoN(-) competition (C and D). Sequencing was performed across 

801 a 1.7kb region of nspl2-RdRp that included both mutations. Data represents 3 independent 

802 experiments, each with 2 replicates. Error bars represent SEM. 

803 
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804 Figure 6. The nspl2-V553I mutation is stable across passage; however, nspl2-M611F is 

805 vulnerable to reversion. DBTs were infected with an initial MOI of 0.01 then blind passaged in 

806 triplicate for 5 passages. Total RNA was taken and sequencing was performed across a 1.7kb 

807 region of nspl2-RdRp that included both mutations. Percentage of each nucleotide present in 

808 each of the triplicate lineages after 5 passages is shown. Mutant viruses in the WT (solid bars) 

809 and nspl4-ExoN(-) (slashed bars) backgrounds are both shown. The original mutation for each 

810 of the viruses is shown above the graph and the likely majority, secondary and tertiary codons 

811 present in the population are shown below the graph. 

812 

813 Figure 7. Resistance of MHV nspl2-RdRp V553I and M611F mutant viruses to 5- 

814 azacytidine in the nspl4-ExoN(-) background. Domain location of mutations is indicated by: 

815 fingers (blues) and palm (reds). DBT cells were pre-treated with different concentrations of 5- 

816 AZC for 30 min. Treatment was removed and cells were infected with indicated viruses at an 

817 MOI of 0.01. Media containing 5-AZC was replaced 30 min p.i. Vims samples were taken at 32 

818 h p.i. and titer was determined by plaque assay. Data represents 5 independent experiments, each 

819 with 2 replicates. Error bars represent SEM (*P < 0.05 by ratio paired t test). 

820 

821 Figure 8. The nspl2-V553I mutation confers decreased accumulation of mutations in the 

822 nspl4-ExoN(-) background with no bias toward the exclusion of specific nucleotides. DBTs 

823 were infected with an MOI of 0.01 and total RNA collected. Deep sequencing was performed on 

824 these samples. The statistically significant mutations present as greater than or equal to 1% of the 

825 total population are shown for wild type, nspl4-ExoN and nspl2-V553I or nspl2-M61 IF in both 

826 backgrounds. These were graphed (A-F) according to their distribution across the genome with 
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827 intentionally introduced mutations (B-F) shown with circles colored blue (nspl2-V553I), red 

828 (nspl2-M61 IF) or green (nspl4-ExoN(-)), as (E) the total number of mutations present in the 

829 population and as (F) the percent of specific mutations present. 
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Table 1. Recovery of mutant viruses using site-directed mutagenesis. 


nspl2-RdRp 

Region 

CVB3 

engineered 
substitutions 
(nspl2) (8) 

MHV 

engineered 

substitutions 

(nspl2) 

Recovery: 

WT (nspl4-ExoN(+)) 

Recovery: 

nspl4-ExoN(-) 

Fingers 

1176V 

V553A 

No 

Not attempted 

V553I 

Yes 

Yes 

Y268H 

Y649H 

Yes 

No 

Y268W 

Y649W 

No 

Not attempted 

Palm 

I230F 

M611F 

Yes 

Yes 

F232Y 

W613Y 

Yes 

No 

A239G 

A621G 

Revertant 

Not attempted 

K360R 

K794R 

Yes 

No 
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